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Abstract — This  paper  considers  the  problem  of  space-time 
adaptive  processing  (STAP)  in  non-homogeneous  environments, 
where  the  disturbance  covariance  matrices  of  the  training  and 
test  signals  are  assumed  random  and  different  with  each  other. 
A  Bayesian  detection  statistic  is  proposed  by  incorporating  the 
randomness  of  the  disturbance  covariance  matrices,  utilizing 
a  priori  knowledge,  and  exploring  the  inherent  Block-Toeplitz 
structure  of  the  spatial-temporal  covariance  matrix.  Speci  cally, 
the  Block-Toeplitz  structure  of  the  covariance  matrix  allows  us  to 
model  the  training  signals  as  a  multichannel  auto-regressive  (AR) 
process  and  hence,  develop  the  Bayesian  parametric  adaptive 
matched  Iter  (B-PAMF)  to  mitigate  the  training  requirement 
and  alleviate  the  computational  complexity.  Simulation  using 
both  simulated  multichannel  AR  data  and  the  challenging 
KASSPER  data  validates  the  effectiveness  of  the  B-PAMF  in 
non-homogeneous  environments. 

Index  Terms — Parametric  adaptive  matched  Iter,  Bayesian 
detection,  space-time  adaptive  signal  processing,  non- 
homogeneous  environments. 

I.  Introduction 

Traditional  space-time  adaptive  processing  (STAP)  usually 
deals  with  homogeneous  environments,  where  the  test  signal  is 
assumed  to  share  the  same  covariance  matrix  with  the  training 
signals  [  1]— [4]. 

To  account  non-homogeneous  environments,  a  number  of 
models  have  been  proposed.  One  is  the  partially  homogenous 
environment,  which  assumes  the  training  signals  share  the  co- 
variance  matrix  with  the  test  signal  up  to  an  unknown  scaling 
factor  [5],  [6],  This  model  can  be  considered  as  a  special  case 
of  the  generalized  eigenrelation  (GER)  [7].  Another  one  is 
the  compound-Gaussian  model,  which  assumes  the  training 
signals  are  a  product  of  a  texture  (scaler)  and  a  Gaussian 
vector.  The  texture  is  used  to  simulate  power  differences 
among  the  signals  from  range  bins  [8],  [9],  More  recently, 
a  new  class  of  non-homogeneous  environments  for  adaptive 
signal  detection  emerges.  This  non-homogeneous  environment 
is  characterized  by  treating  disturbance  covariance  matrices  of 
both  the  test  signal  and  training  signals  as  random  matrices 
and  ensuring  that  they  are  different  in  probability  one  [10]— 
[14].  Following  the  Bayesian  non-homogeneous  model  in  [10], 
an  adaptive  matched  Iter  has  been  derived  by  replacing  the 
exact  covariance  matrix  of  the  test  signal  by  its  maximum  a 
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posteriori  (MAP)  estimate  in  the  matched  Iter  [10].  It  has 
been  shown  that,  by  accounting  the  heterogeneity  knowledge, 
the  Bayesian  adaptive  matched  Iter  (B-AMF)  outperforms 
the  standard  AMF  in  the  non-homogeneous  environment.  For 
applications  of  this  Bayesian  non-homogeneous  model  that 
employ  space-time  adaptive  processing,  the  training  require¬ 
ment  of  the  sample  covariance  matrix  (SCM)-based  B-AMF 
cannot  be  met  due  to  the  scarcity  of  the  training  signal.  For 
example,  with  J  =  11  spatial  channels  and  N  =  32  coherent 
pulses  of  the  KASSPER  dataset  and  assuming  an  instantaneous 
RF  bandwidth  of  500  KHz,  K  =  JN  =  352  training  signals 
calls  for  the  training  range  over  a  200-km  range,  which  is  not 
practical  [15],  [16], 

In  this  paper,  while  preserving  the  Bayesian  non- 
homogeneous  environment,  we  further  explore  the  inherent 
Block-Toeplitz  structure  of  the  spatial-temporal  covariance 
matrix  which  allows  the  block  LDU  decomposition  [17],  and 
hence  enables  the  disturbances  to  be  modeled  as  a  multi¬ 
channel  auto-regressive  (AR)  process  [17]-[23].  The  result¬ 
ing  Bayesian  parametric  adaptive  matched  Iter  (B-PAMF) 
reduces  the  joint  spatial-temporal  whitening  of  the  SCM- 
based  B-AMF  to  successive  spatial  and  temporal  whitening. 
As  a  result,  it  facilitates  the  STAP  in  the  non-homogeneous 
environments,  and  reduces  the  excessive  training  requirement 
of  the  B-AMF.  Moreover,  the  B-PAMF  is  able  to  incorporate 
heterogeneities  of  the  signals  and  utilize  available  a  prior 
knowledge  to  the  decision  statistic.  The  effectiveness  of  the 
B-PAMF  is  veri  ed  by  using  the  simulated  multichannel  AR 
data  and  the  high  delity  KASSPER  data  [15]. 

II.  Signal  Model 

Assume  J  spatial  channels,  N  temporal  pulses,  and  K 
training  range  cells.  The  problem  of  interest  is  to  detect  a 
JN  x  1  multichannel  signal  s  with  unknown  amplitude  a  in 
the  presence  of  spatially  and  temporally  correlated  disturbance 

d0: 

H0-  x0(n)  =  d0(n),n  =  0, 1,  •  •  •  ,  N  -  1, 

Hi  :  xo (n)  =  as(n )  +  do(n),  n  =  0, 1,  •  •  •  ,N  —  1, 

where  do  =  [dj (0),  dj (1),  •  •  •  ,  dj ( N  —  1)]T,  and  s  and  xo 
are  similarly  de  ned.  In  this  paper,  the  signal  model  makes 
the  following  assumptions: 

•  AS1  (Multichannel  AR  Process):  The  disturbances  in 
both  test  and  training  signals  are  modeled  as  a  multi¬ 
channel  AR  process  [17] — [19]: 

dfc(n)  =  -  J2i= i  AH(i)dk(n-i)+ek(n),  k  =  0,  •  •  •  ,K, 

(2) 
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where  AH  =  [Ap(l),  Ap (2),  •  •  •  ,  A H(P)]  denote  the 
unknown  multichannel  AR  coef  cient  matrix,  and  £k(n) 
denote  the  J  x  1  temporally  white  but  spatially  colored 
noise  vectors. 

•  AS2  (Random  Disturbance  Covariance  Matrix  of 
Training  Signals):  The  noise  vector  £fc(n)  is  distributed 
as  £k{n)  ~  CJ\f{ 0,  Q),  and  the  spatial  covariance  matrix 
Q  follows  an  inverse  complex  Wishart  distribution  with 
degrees  of  freedom  p  and  mean  Q: 


p(  Q) 


J)Qr  e-(M-J)  tr(Q_1Q) 

f(J, p)  |Q|(/i+J) 


(3) 


where  f  ( J,  /r)  =  7rJ(J_1)/2  11^=1  r  {p  —  J  +  fc)  with  T 
given  by  the  Gamma  function  [24]. 

•  AS3  (Random  Disturbance  Covariance  Matrix  of  Test 
Signal):  The  noise  vector  in  the  test  signal  £o(n)  ~ 
CAT( 0,  Qo),  and  Qo,  given  Q,  has  a  complex  Wishart 
distribution  with  degrees  of  freedom  v  and  mean  Q 


p(Qo|Q) 


^IQor- V  ^trfQ-tQo) 
f  (j,  v)\Q\v 


(4) 


The  multichannel  AR  process  for  the  disturbances  in  both 
the  test  and  training  signals  consists  of  two  types  of  unknown 
parameters:  one  is  the  deterministic  AR  coef  cient  matrix  A, 
and  the  other  is  the  random  spatial  covariance  matrices  Q 
and  Qo.  The  available  a  prior  knowledge  is  imposed  on  the 
mean  of  Q,  i.e.,  Q,  which  can  be  obtained  from  sources 
such  as  land-use  maps,  past  measurements,  etc  [15].  It  is 
said  that  the  importance  of  the  a  priori  knowledge  Q  is 
controlled  by  parameter  //,  while  the  the  heterogeneities,  i.e., 
the  statistical  differences  between  the  test  and  training  signals, 
are  determined  by  parameter  v.  Most  importantly,  Q  /  Q0 
with  probability  one,  which  ensures  the  environment  non- 
homogeneous  [10]. 


III.  Bayesian  Parametric  Adaptive  Matched  Filter 

By  assuming  perfect  knowledge  on  parameters  A  and 
Qo,  the  solution  to  the  problem  of  interest  is  the  classical 
parametric  matched  Iter  (PMF)  [17]: 

N—l  2 

E  §ff(n)Qo1x0(?i)  Ri 

Tpmf  =  n~NP_1 -  7pmf,  (5) 

E  sH(n)Qo  1s(n)  H° 

n—P 


statistic  by  their  estimates  (e.g.,  the  maximum  likelihood 
estimate  (MLE)  by  using  the  training  signals).  For  the  non- 
homogeneous  environment  considered  in  this  paper,  due  to  the 
randomness  of  the  spatial  covariance  matrix,  we  adapt  a  hybrid 
parametric  AMF,  which  is  denoted  as  the  Bayesian  PAMF  (B- 
PAMF)  by  rst  obtaining  the  MLE  of  the  deterministic  AR 
coef  cient  matrix  A  ml,  then  deriving  a  maximum  a  posteriori 
estimate  (MAP)  of  the  stochastic  spatial  covariance  matrix 
Qo.map  and  nally  replacing  A  and  Qo  in  the  PMF  statistic 
by  their  estimates. 


A.  MLE  of  A 

According  to  the  signal  model,  the  joint  probability  density 
function  (pdf)  of  the  training  signals  can  be  approximated 
(ignore  the  conditionality  on  the  rst  P  temporal  vectors)  as 
[18] 


/(x i,---  ,xa-|A,Q)  = 
where 


1  r-tr(Q-1S(A)) 

7TJ  |Q  j 


1  K(N-P) 


K  N—l 


S(A)  =  K(N  -  P )  E  E  £k(n)e%(n) 


(8) 


k—1 n—P 


From  AS2,  we  can  remove  the  dependence  of  the  above  pdf 
on  Q  by  integrating  it  over  Q: 


/(x  !,-•*  ,xa|A)  =  y  /(x  1, 

\(P-mr  p 


,  xa|A,  Q)p(Q)dQ 


nJK(N-P)Y ( j  ^ 

| (ft-  J)Q|A1f(J,L-  J) 

nJK(N-P)f  ( J  ^ 


QrLe-tr(Q"lf:)dQ 

S|'i+J  (9) 


where  L  =  p  +  J  +  K(N  -  P)  and  S  =  K{N  -  P)S(A)  + 
(p  —  J)  Q.  Therefore,  nding  the  MLE  of  A  is  equivalent  to 
minimizing  the  determinant  of  S.  Rewrite  the  matrix  S  as 


S  =  K(N  -  P)S(A)  +  (p-  J)Q 
=  R-a-a:  +  R-yx  A  +  AHjk,yx  +  A^Ry^A  +  (p  —  J)  Q 

=  Ryy  ^A^  +  R^Ryy^ 

+  R-**  —  +  (fl  —  j)  Q)  (10) 


where 


where  7pmf  denotes  the  PMF  threshold  subject  to  a  selected 
probability  of  false  alarm,  and  the  whitened  steering  vector 
and  test  signal  are  obtained  by  using  the  true  AR  coef  cient 
matrix  A 

p 

s(n)  =  s(n)  4-  ^  AH (p)s(n  -  p),  (6) 

p=  i 
p 

x0(n)  =  x0()t)  +  ^  AH(p)x0(n  -  p).  (7) 

P=  1 

The  parametric  AMF  replaces  the  exact  AR  coef  cient  ma¬ 
trices  A  and  the  spatial  covariance  matrix  Qo  in  the  PMF 


K  N- 1  K  N- 1 

=  Xfc(n)xf  (n),Ryy  =  E  y  fc(n)yf(n)> 

k= 1 n=P  k—1 n=P 

K  N-l 

^^EE^w^W’  (11) 

k—1 n—P 

with  y k(n)  =  [x^(n  -  1),  •  •  •  ,  x^ (n  -  P)] T  G  CJPxl. 

Since  R^y  is  nonnegative  de  nite  and  the  remaining  terms 
R.IX  —  Rp  RyyRyX  +  (p  —  J)Q  do  not  depend  on  A,  it 
follows  from  (10)  that 

S  ^  S|  A=AMl  =  ftra:  ~  R-yx^-yy  ^yx  +  (f*  —  ^)Qj  (12) 
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where  the  MLE  of  A  is  given  as 

AmL  =  —  R^Ryy-  (13) 

When  S  is  minimized,  the  MLE  Aml  will  minimize  any 
nondecreasing  function  including  the  determinant  of  S. 


B.  MAP  Estimate  of  Qo 

The  MAP  estimate  of  Qo  requires  the  computation  of  the 
posterior  distribution  /  (Qo|xi,  X2,  •  •  •  ,  x^): 

/  (Qo|xi,x2,  •  •  •  ,x/f)  =  /  /  (Qo>  Q|xi,x2,  •  •  •  ,xK)dQ, 

(14) 

where 


/(Qo,Q|xi,x2,--.:  ,xk) 

oc/(xi,x2,>  •  •  ,x*r|Q0,  Q)p(Qo|Q)p(Q) 

a  |Qor_J  |Q|-(i+y)  e-te(Q“1[=+«'Qo]). 


As  a  result,  (14)  can  be  calculated  as 
J f  (Qo>  Qlxi,  x2,  •  •  •  ,xK)dQ 

<x  |Qor_J  J  |Qr(i+!/) 

a  |Qor_J  |S  +  ^Qo|M+I/+^(JV“jP).  (16) 

Taking  the  logarithm  of  the  above  equation,  then  taking  the 
derivative  with  respect  to  Qo,  and  equaling  to  zero,  we  have 

dln/(Q0|xi,x2,  •  •  ■  ,xK) 
dQo 

oc  [v  -  J)Qq  1  -  v  (p  +  v  +  K(N  -  P))  [S  +  ^Qo]^1  =  0, 

(17) 


which  suggests  that,  given  A,  the  estimate  of  Qo  is 

(v-J) 


Qo  — 


v(n+J+K(N-P)) 


£. 


(18) 


where  7b-pamf  denotes  the  B-PAMF  threshold  subject  to  a 
selected  probability  of  false  alarm,  Qo.map  is  given  by  (19), 
and  the  whitened  steering  vector  and  the  whitened  test  signal 
are  obtained  by  using  Aml  given  by  (13): 

p 

s (n)  =  s (n)  +  AmlO)s(>i  ~p),  (21) 

P=  1 
P 

x0(n)  =  x0(n)  +  Aml(p)xo(«  -  p)-  (22) 

p=  i 

From  (20),  on  one  hand,  it  is  seen  that  the  B-PAMF  per¬ 
forms  successive  whitening,  i.e.,  temporal  whitening  followed 
by  spatial  whitening,  as  opposed  to  joint  spatio-temporal 
whitening  across  all  JN  dimensions  of  the  Bayesian  AMF 
in  [10],  On  the  other  hand,  compared  with  the  standard 
PAMF  [17],  the  B-PAMF  incorporates  the  a  priori  knowledge, 
i.e.,  Q,  and  utilizes  the  heterogeneity  parameter  v  and  the 
importance  parameter  of  the  a  priori  knowledge,  i.e.,  p  into 
the  estimate  of  the  spatial  covariance  matrix  Q.  Hence,  it  is 
allowed  that  the  B-PAMF  provides  computational  ef  ciency 
and  mitigates  training  requirement,  meanwhile  exploiting  the 
a  priori  knowledge  and  the  heterogeneity  to  improve  the 
performance  of  detection. 

IV.  Numerical  Examples 

In  this  section,  simulation  results  are  provided  to  illustrate 
the  performance  of  the  B-PAMF.  Speci  cally,  we  rst  test 
the  B-PAMF  by  using  simulated  data  which  conforms  to 
AS1,  AS2  and  AS3,  and  then  using  the  more  challenging 
KASSPER  dataset  [15].  The  disturbance  signal  is  generated 
as  a  multichannel  AR(2)  process  with  AR  coef  cient  A  and  a 
spatial  covariance  matrix  Q.  The  signal  vector  s  corresponds 
to  a  uniform  equispaced  linear  array  with  randomly  selected 
normalized  spatial  and  Doppler  frequencies.  The  signal-to- 
interference  plus  noise  ratio  (SINR)  is  de  ned  as 

SINR  =  lofs^R^s,  (23) 


Replacing  A  with  Aml  of  (13)  in  the  above  estimate  (viz,  S), 
the  MAP  estimate  of  Qo  is 


Qo.map  — 


(u- J) 


v(n  +  J  +  K{N  -  P)) 

R*z  —  Ryz-RyyRj/z  +  (p  ~  <7)Q 


(19) 


It  is  seen  that  the  MAP  of  Q  is  a  linear  combination  of  a 
standard  estimate  of  Q  as  introduced  in  [18],  [19]  and  the  a 
priori  knowledge  Q.  This  linear  combination  has  been  seen 
before  for  non-parametric  approaches  [10],  [25]. 


C.  Bayesian  PAMF 

With  the  ML  estimate  of  A  and  the  MAP  estimate  of 
Qo,  the  adaptive  version  of  the  PMF  in  the  heterogeneous 
environment  can  be  derived  as 

N- 1,  „  „  2 

E  sff(n)Qo,MAP*o(n)  Ri 

7b-PAMF  =  x -  s;  7B-PAMF  (20) 

E  §ff(n)Qo,MAPi(rl)  H° 

n—P 


where  R  corresponds  to  the  assigned  AR  coef  cient  matrix  A 
and  the  mean  of  the  training  spatial  covariance  matrix  Q.  For 
each  Monte-Carlo  trial,  the  spatial  covariance  matrix  Q  for  the 
training  signal  is  generated  as  an  inverse  Wishart  distribution 
with  mean  Q,  and  then,  given  Q,  the  spatial  covariance  matrix 
Qo  for  the  test  signal  is  generated  as  a  Wishart  distribution 
with  mean  Q. 

We  focus  here  on  performance  comparison  between  the  B- 
PAMF  and  the  standard  PAMF  [17]  in  the  non-homogeneous 
environment.  Figs.  1(a)- 1(c)  show  the  probability  of  detection 
versus  the  SINR  for  the  B-PAMF  and  the  standard  PAMF  in 
cases  of  different  values  of  /j  and  v,  when  P  =  2,  J  =  4, 
N  =  16,  K  =  1,  and  Pf  =  0.01.  It  is  seen  that,  in  all 
three  cases,  the  B-PAMF  outperforms  the  standard  PAMF. 
Speci  cally,  for  a  xed  value  of  p,  a  larger  value  of  /y, 
more  homogeneous  environments,  results  in  slightly  improved 
performance  of  detection  of  the  B-PAMF,  from  Fig.  1(a)  to 
Fig.  1(b).  On  the  other  hand,  for  a  xed  value  of  v,  increasing 
p  means  more  importance  of  the  a  priori  knowledge  of  Q, 
which  leads  to  wider  performance  gap  between  the  B-PAMF 
and  PAMF,  as  compared  between  Fig.  1(a)  and  Fig.  1(c). 
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K=1 ,  J=4,  N=16,  P  =0.01,  AR(2),  jx=9,  v=5  K=1,  J=4,  N=16,  P  =0.01,  AR(2),  ^=9,  v=10  K=1,  J=4,  N=16,  P=0.01 ,  AR(2),  ^=14,  v=5 


(a)  (b)  (c) 


Fig.  1.  Probability  of  detection  versus  SINR  (a)  when  P  =  2,  p  =  9,  and  v  =  5;  (b)  when  P  =  2,  p  =  9,  and  v  =  10;  (c)  when  P  =  2,  p  =  14,  and 
i/  =  5. 


V.  Conclusion 

A  Bayesian  parametric  adaptive  matched  Iter  has  been 
proposed  by  modeling  the  disturbances  in  the  test  and  training 
signals  as  a  multichannel  AR  process  and  simulating  the  het¬ 
erogeneity  between  the  training  and  test  signals  by  introducing 
random  disturbance  covariance  matrices.  The  B-PAMF  admits 
successive  temporal  and  spatial  whitening,  which  reduces 
the  computational  complexity  of  the  joint  spatial-temporal 
whitening  based  adaptive  detectors.  The  training  requirement 
is  also  reduced.  Simulation  results  validate  the  effectiveness 
of  the  B-PAMF. 
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